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Abstract 

Azimuthal spin asymmetries, both for charged and neutral pion production in semi- 
inclusive deep inelastic scattering of unpolarized charged lepton beams on longitudi- 
nally and transversely polarized nucleon targets, are analyzed and calculated. Various 
assumptions and approximations in the quark distributions and fragmentation func- 
tions often used in these calculations are studied in detail. It is found that different 
approaches to the distribution and fragmentation functions may lead to quite differ- 
ent predictions on the azimuthal asymmetries measured in the HERMES experiments, 
thus their effects should be taken into account before using the available data as a 
measurement of quark transversity distributions. It is also found that the unfavored 
quark to pion fragmentation functions must be taken into account for tt~ production 
from a proton target, although they can be neglected for 7r+ and 7r° production. Pion 
production from a proton target is suitable to study the u quark transversity dis- 
tribution, whereas a combination of pion production from both proton and neutron 
targets can measure the flavor structure of quark transversity distributions. 



PACS numbers: 13.87.Fh, 13.60.-r, 13.88. +e, 14.20.Dh 



1 Introduction 



Recently, the HERMES collaboration reported the observation of single-spin az- 
imuthal asymmetries for charged and neutral pion production, in semi-inclusive deep 
inelastic scattering (DIS) of unpolarized positron beam on the longitudinally polar- 
ized proton target Q. Such azimuthal asymmetries are important, because they 
can provide information of the transversity distributions , which are one of the three 
fundamental quark distributions of the nucleon. The other two fundamental quark 
distributions are the momentum and helicity distributions, which have been explored 
for more than three decades in various DIS processes @, ^, and their explicit flavor- 
dependence has been obtained with relative high precision, although there are still 
some uncertainties concerning the sea content at small Bjorken variable x and the va- 
lence flavor-spin structure as x 1. However, the experimental measurements of the 
quark transversity distributions are just about to begin, since the transversity is not 
directly observable in inclusive DIS processes. It has been proposed that the transver- 
sity can manifest itself through the Collins effect ||^ of nonzero production between a 
chiral-odd structure function and a T-odd fragmentation function, which is accessible 
in some specific semi-inclusive hadron production experiments [^, |^, H, |10|, 11, |l^ . 
Indeed, there have been a number of studies [0, [T|, 0, |TB|, [T^, 0, ig ^, 0, 0, ^, ^ 
to show that the azimuthal asymmetries measured by HERMES can be used to ex- 
tract the quark transversity distributions of the nucleon, although the theoretical 
predictions of the explicit flavor-dependent transversity distributions are different in 
various models |^ . Therefore the HERMES measurements provide us the opportunity 
to confront theoretical predictions with physical observations, and it is meaningful to 
examine the theoretical results carefully, in order to see what conclusion we can draw 
on the transversity. 

The analyzing power of azimuthal spin asymmetry measured by HERMES is de- 
fined as 

,w ^ /[d0]iy(0) {iv+(0)-iv-(0)} 

i/[d0]{iV+(0) + iV-(0)} ' ^ ^ 

where UL denotes unpolarizedheam on a longitudinally polaiized target, W{(j)) = sm(j) 
or sin 20 is the weighting function for picking up the Collins effect, and N~^{(f)) {N~{(j))) 
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is the number of events for pion production, as a function of 0, when the target 
is positively (negatively) polarized. The azimuthal angle is the angle between 
the pion emitting plane and the lepton scattering plane, with the lepton scattering 
plane determined by the incident and scattered leptons, and the pion emitting plane 
determined by the final detected pion and the virtual photon. The virtual photon acts 
as the common axis of both planes. The analyzing powers of azimuthal asymmetries 
for charged pions, tt~^ and vr", and neutral pion, vr", have been measured [|l], 0], and 
there is clear evidence for non-zero values of A^^l^^ for 7r+ and tt'^ production, which 
indicates the existence of azimuthal asymmetries. Under a number of simplifying 
assumptions and approximations, the analyzing power A^^lf' of azimuthal asymmetry 
is proportional to the ratio 

T.,elq{x)Dl{z) ' ^ ' 

Here Cg is the charge of the quark with flavor g, q{x) and 5q{x) are the quark momen- 
tum and transversity distributions of the nucleon target, D\{z) is the fragmentation 
function for an unpolarized quark with flavor a into a pion tt, and Hi^^'^'^{z) is the 
so called Collins function describing the fragmentation of a transversely polarized 
quark into a pion tt. If one assumes, as adopted in almost all previous analysis, that 
the pion production is dominated by the "favored" fragmentation processes u tt^ 
and c? — s> TT^ for charged pions, and neglects the "unfavored" processes u ^ and 
(i 7r+, one may conclude that the 7r+ azimuthal asymmetry measures the -u-quark 
transversity 6u{x) and the 7i~ azimuthal asymmetry measures the d-quark transver- 
sity 6d{x), respectively. This is the situation of the current available studies. 

The purpose of this paper is to take into account the terms neglected in (0), and 
also to take into consideration the contribution from the "unfavored" fragmentation 
process. We will estimate their effects on the azimuthal asymmetries for the three 
kinds of produced pions: 7r+, tt^, and vr''. We find that the neglected terms may have 
sizable effects on the azimuthal asymmetries in the HERMES kinematical region, and 
that the results are sensitive to different assumptions on the quark distributions and 
fragmentation functions. We find also that, although the neglect of the "unfavored" 
fragmentation processes seems to be reasonable for 7r+ production, it is not justi- 
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fied for tt production. Moreover, a significant difference between the vr azimuthal 



asymmetries predicted by two different models |^0[; ^ire largely reduced by taking into 
account the "unfavored" processes, and this renders it difficult to use tt^ production 
as a clean process to test different predictions on the flavor-spin quark structure of 
the nucleon at large x. We will also check the possible ambiguities from the uncer- 
tainties in the "unfavored" fragmentation functions. However, it will be shown that 
the favored fragmentation process in Tr~^ production from a neutron target plays a 
more important role than that in the corresponding tt^ production from a proton 
target. Thus tt^ production from a neutron target can provide a test of different 
model predictions. A combination of pion production from both proton and neutron 
targets can measure both of the u and d quark transversity distributions, 6u{x) and 
6d{x), of the proton. 

This paper is organized as follows. In Section II of the paper, we will present the 
formulae for the azimuthal asymmetries, and analyze in details various assumptions 
and approximations on the quark distributions and fragmentation functions used 
for numerical calculations. In section III, we perform numerical calculations of the 
azimuthal asymmetries for pions: 7r+, tt~ , and vr*^, and give detailed discussions on the 
effects of using different approaches to the distribution functions and fragmentation 
functions. We discuss also in detail the uncertainties coming from the "unfavored" 
processes of quark to pion fragmentation. In Section IV, we predict the azimuthal 
asymmetries for pion production in semi-inclusive DIS processes of an unpolarized 
charged lepton beam on a transversely polarized nucleon target. Finally, we present 
the summary and conclusions in Section V. 



2 Formalism of Azimuthal Spin Asymmetry in Pion 
Electroproduction 

Deep inelastic scattering (DIS) of a lepton on a nucleon is one of the most efficient 
and clean tools to detect the underlying structure of the nucleon. With the recent 
progress in experimental techniques, the precision detection of hadrons produced in 
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DIS processes makes it possible to reveal more about the basic structure of the nu- 
cleon, and the use of polarized beams and/or polarized targets makes it also possible 
to measure with some precision various polarized quark distributions and fragmenta- 
tion functions. The HERMES experiments [jl], ^ are performed with an unpolarized 
lepton beam scattering on a longitudinally polarized proton target, and some assump- 
tions and approximations are needed in order to connect the pion production with 
the transversity distribution, which is supposed to be related to the probability of 
finding a quark with its spin aligned to the proton spin for a transversity polarized 
proton. We first give the formulae for the analyzing power of azimuthal asymmetry 
and analyze the quantities needed to perform numerical calculations. 



2.1 Spin asymmetries of pion electroproduction 

We now introduce the basic kinematics for semi-inclusive DIS processes. Let h and 
I2 be the initial and final momenta of the incoming and outgoing charged leptons 
respectively, = —q^ be the squared 4-momentum transfer of the virtual photon 
with momentum q = I1 — I2, and P and Ph (M and M^) be the target and final hadron 
momenta (masses). Then the three basic variables, x, y, and z, used in describing 
DIS processes can be expressed as 

Q' P-q P-Ph 

2P ■ q P ■ k P-q 

where x is the Bjorken variable which corresponds to the momentum fraction carried 
by the struck quark in a light-cone description, y is the fraction of the initial lepton's 
energy transferred to the quark by the virtual photon, and z is the fraction of quark 
momentum transferred to the final produced hadron in the fragmentation process. In 
the target rest frame, the three variables can be expressed as 

where z/, and Eh are the virtual photon energy, the incident lepton energy, and the 
produced hadron energy, respectively, s = 2P-1 = 2ME is the total squared energy in 
the lepton-proton center of mass frame, therefore we have a relation Q"^ = sxy. The 
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incident lepton energy is = 27.6 GeV in the HERMES experiments. The invariant 
mass of the photon-proton system is 1^ = \/1Mv + — Q^. Let Phi. be the final 
hadron transverse momentum and l\i_ be the initial lepton transverse momentum, 
both respect to the virtual photon momentum direction. Then the azimuthal angle 
is the angle between P^x. and around the virtual photon direction. 

Using the above measurable variables of semi-inclusive hadron production and the 
measurable cross sections, one can defined sin and sin 20 weighted spin asymmetries 

' Sm0 = r^9n <^.^ , ^-N ' (5) 



\ Mn "/ /d2p^^(da+ + da-) 

\PkA^ . , ,\ /d^-P/^T§tfsin20(da+-da-) 



\MMh ^1 jd^Ph^{da+ + da- 

which are measurable quantities. Here -|-(-) denotes parallel (antiparallel) longitudi- 
nal polarization of the target. The above asymmetries are related to the azimuthal 
asymmetries measured by HERMES through the following relations 

A sin 4, ^ 2M^ l\Phi_\ . \ 

' sin 20). (8) 



h±\ 



i,2</, _ 2MMfe / |P; 

~ {PL) \MM, 



In the HERMES experiments, {Ph±) = 0.44 GeV, and we may use (P^j_) = (0.44)^ 
GeV^ as an estimate of {Ph±) in our calculations. 

Prom another side, one can use the quark-parton model to calculate hadron pro- 
duction in semi-inclusive DIS processes, and express the cross sections in terms of 
various distribution and fragmentation functions. Then one can define the weighted 
cross section for one fiavor quark, dag, with different weighting functions depending 
on the final hadron transverse momenta Wi{Ph±): 

Ef = J d^Ph±Wi{Ph±)dag, (9) 

where Wi{Ph±) = 1, Wi{Ph±) = Ph±sm(P/Mh, and w,{Ph±) = \Ph±\^ sm24>/MMh. 
Thus one can relate the theoretical calculations with the above spin asymmetries by 



where a sum over all quark flavors, Sj = X^nC^S^, is implicitly assumed, and will be 



assumed from now on. In the case of unpolarized beam and longitudinal polarized 
target. Si, S2 = S2L + S2T, and S3 are given by |0, ||, |2 



Si= l + {l-yf fi{x)D,{z), 



M 



Q 



S2L = 4^L^(2 - y)Jl - y xhL{x)zH^^'\z) - htl'\x)H{z 



S2T = 2St,{1 - y)h^{x)zHt^^\z) 



Jls = 8Sl{1 - y)hti'\x)z'H^-[z). 
Here the components of the longitudinal and transverse target polarization in the 



-(1), 



(12) 

(13) 

(14) 
(15) 



virtual photon frame are denoted by Sl and Stx, respectively: 



Sl = Scos9^, Stx = •S' sin ft 



7' 



(16) 



with 



sin 6^ 



1-y- 



X 



M^xy' 



\ sy + AM^ 



X 



(17) 



and S is the target polarization, which has the value 0.86 in the HERMES experi- 
ments. Twist-2 distribution functions and fragmentation functions have a subscribe 
"1": /i and Di are the usual unpolarized distribution and fragmentation function, 
while hf[^\x) and hi{x) are the quark transverse spin distribution functions of lon- 
gitudinally and transversely polarized nucleons, respectively. /il(x) is the twist-3 
distribution function of a longitudinally polarized nucleon, and it can be split into a 
twist-2 part, hfj^^\x), and an interaction dependent part, hiix): 



hiix) 



h 



-(1), 

IL ' 



X] 



X 



+ hL{x). 



The fragmentation function H{z) is the interaction dependent part of the twist-3 
fragmentation function: 



H{z) = -2zHt^^\z)+H{z). 



(19) 



The functions with superscript "(1)" denote p^- and /c^-moments, respectively: 



Hti'\z) ^ z'Jd^k^^^H^,{z,z^kl), (21) 

where p± and k± are the intrinsic transverse momenta of the initial and final partons 
in the target and produced hadrons, respectively. 

To calculate the spin asymmetries and compare them with experiments, we need 
to know the quark distribution functions: fi{x), hi{x), hiix), and hf[^\x), and the 
fragmentation functions: Di{z), Hi^^\z), and H{z). Unfortunately, most of the 
distribution functions and fragmentation functions which appear in these expressions 
are not known a priori, since they have not been measured yet. Thus we have to 
make some assumptions and approximations, and this leads to different approaches 
for the distribution functions and fragmentation functions [|T^ 



2.2 Different approaches to the distribution functions and 
fragmentation functions 

The most simple approach, denoted as Leading Approach, is to neglect the 1/Q term 
Ti2L in ^2, i-e., we neglect both the hiix) and hij^\x) terms in the spin asymmetry 
^^^sin0^. Then we find immediately that 

I\Pki] . A S2T h,{x)Ht'-'\z) 

■sm0) = oc — — — — — , (22) 



\ Si h{x)D,{z) 

which is Eq. (|]). Thus we can relate the quark transversity distribution hi{x) to 
the spin asymmetry with its x-dependence. This is the situation discussed in several 
MM- 



papers 



The next approach, denoted as Approach 1, is to assume that the twist-2 quark 
transverse spin distribution function of longitudinally polarized nucleon, hfl^^\x), is 
zero [^. Then it follows that 



h^x) = hix) = h{x). (23) 
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Another approach, denoted as Approach 2, is to assume that the interaction de- 
pendent twist-3 part, hi^x), is zero and consequently, we can also assume that 



H{z) is zero. Then by neglecting the term proportional to the current quark mass, 
one can obtain a Wandzura-Wilczek type relation |2^, |26| 



fti-f(x) = -x=/'dfM), (24) 



X 



hL{x) = 2xl de^. (25) 



It follows, from Eq. (|T|), that 

^1 Mi) 

e 

In the above two approaches, we only need the distribution functions, and 
hi{x), and the fragmentation functions, Di{z) and Hi^'^\z), in order to calculate the 
spin asymmetries (^^^sin0^ and ^^^^sin! 



2.3 Quark transversity distributions of the nucleon 

In the three approaches described above, we need the quark distribution functions, 
= q{x) and hi{x) = Sq{x), as inputs. The x-dependence of the spin asymmetries 
are controlled by the above two quark distribution functions, in which the unpolarized 
quark distributions q{x) are known to relative high precision, thus we are able to 
extract the information on the quark transversity distributions 6q{x) by comparing 
theoretical predictions with the experimental data on spin asymmetries. 

Although the quark momentum distributions q{x) have been measured to high pre- 
cision, there are still uncertainties concerning the flavor structure as a; ^ 1. For exam- 
ple, the SU(6) quark spectator diquark model [|^, ^ ^ predicts d{x)/u{x)\,^^-^ = 0, 



whereas a pQCD based counting rule analysis |^, |3T|, ^] predicts d{x)/u{x 



\x=l 

1/5. The flavor structure of the transversity distributions in the two models are also 
found to differ significantly in this regionpO, |33[: the pQCD based analysis predicts 



20| , |33| 6d{x)/d{x) 1, whereas the SU(6) quark-spectator-diquark model predicts 



3J, ^ 6d{x) / d(x) —1/3. Detailed constructions of the quark transversity distri- 



butions in the two models can be found in Refs. |33 |. 

To make realistic predictions of measurable quantities, we need also to take into 
account the sea quark contribution in the two model constructions. In the quark 
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diquark model case, this can be achieved by adopting one set of unpolarized quark 
distribution parametrization as input, and then use theoretical relations to connect 



the quark transversity distributions with the unpolarized quark distributions ||29|, |4 



5u^{x) = [u^{x) - ld^{x)]Wsix) - ld^{x)Wvix); ^^^^ 
6dv{x) = —\dy{x)Wv{x), 



in a similar way as was done for the quark helicity distributions p9|. Ws{x) and 



Ws{x), which come from the relativistic effect of quark transversal motions are 



the Melosh-Wigner rotation factors |35|, qT^ for spectator scalar and vector di 



quarks. We use the valence quark momentum distributions u^{x) and dy{x) from 
quark distribution parametrization, but with Ws{x) and Wv{x) from the model cal- 
culation [Q. In this way we can take into account the sea contribution for the 
unpolarized quark distributions through the input parametrization. This can provide 
a more reliable prediction for the absolute magnitude and shape of a physical quantity 
than directly from the model calculation. 

For the pQCD based analysis, we adopt the same consideration as above and 
make the following connection to relate the pQCD model quark distributions with 
the parametrization 



nf''''{x) = ur%x), df^'^ix) = %f(«r™(x), (27) 



Suf^^ix) = ^^«r"(x), 6df^-^ix) = ^|/^«r"(^), (28) 

"U^ [X) [^X) 

where the superscripts "th" means the pure theoretical calculation in the pQCD 
analysis |3^ , and "para" means the input from parametrization. The superscript 



"pQCD" means that the new quark distributions keep exactly the same flavor and spin 
structure as that in the pQCD analysis, but with detailed x-dependent behaviors more 
closed to the realistic situation. This is equivalent to use an unique correction factor, 
^para^^y^th^ to adjust cach pure theoretically calculated quantity to a more realistic 
pQCD model quantity. In this way we can take into account the sea contribution by 
using the sea quark distributions from parametrization, while still keep the pQCD 
model behaviors of the valence quark distributions. 
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Figure 1: The different ffavor and transversity structure of the two model quark 
distributions: (a) dv{x)/uv{x)] (b) 5uy{x) /u^{x)] and (c) dd^i^x) / d^{x). The sohd 
and dashed curves correspond to the quark-diquark model and the pQCD analysis, 
respectively. 

The explicit flavor and transversity structure of the two sets of quark distributions 
can be found in Fig. |l]. Thus we have two sets of quark distributions of q{x) and 5q{x), 
which keeps the same valence behaviors as in the quark diquark model and pQCD 
based model predictions, but with the sea from parametrization. We expect to make 
some more realistic predictions based on these two sets of quark distributions, while 
still can distinguish between the different flavor and spin structure predicted by the 
two models. 

There might be some uncertainties of the magnitude of the explicit quark dis- 
tributions, however, these uncertainties will be reduced since we will only calculate 
quantities which involve the ratios of different quark distributions in our later numer- 
ical calculations. 



2.4 Fragmentation functions of pion 

We now turn our attention to the fragmentation functions Di{z) and Hi^^\z). Di{z) 
is the unpolarized fragmentation function for a quark splitting into a pion, and in pre- 
vious studies on the azimuthal asymmetries, only the favored fragmentation functions, 
e.g., for a, u oi d quark into tt^, are considered. However, a new parametrization 
of quark to pion fragmentation functions has been proposed [Q, which combines 
pion production data from both semi-inclusive DIS and e^e~ processes. This new 
parametrization provides a complete set of both favored and unfavored fragmenta- 



11 



tion functions, and therefore it is very useful for the extraction of nucleon structure 
information from pion production in semi-inclusive DIS processes. In this paper, we 
will adopt this set of quark to pion fragmentation functions in our analysis of the 
transversity distributions through the azimuthal spin asymmetries discussed above. 

For the favored fragmentation function, we mean the fragmentation function de- 
scribing the transition of a "valence" quark into the pion, and it is, assuming isospin 
symmetry, 

D{z) = Df{z) = Df{z) = Dliz) = Df{z) (29) 

for charged pion fragmentation. The unfavored fragmentation function is supposed 
to describe the transition of a light-flavor "sea" quark into pions, and it is 



D{z) = Dl\z) = Df{z) = Df{z) = Dl{z) 



(30) 



for vr fragmentation. For the case of the neutral pion, assuming that 



d; (z) + d; (z) 



we find 



(31) 



D-\z) = D:\z) = Dl\z) = Df{z) = Dfiz) 



D{z) + b{z) 



(32) 



The explicit analytically forms of D{z) and D{z) are given by |]38 



/^(Z) = 0.689^-1-039(1 -2)1-241^ 
D(;z) =0.217^-1-805(1 -2)2.037_ 



(33) 



From another point of view, assuming the Gribov-Lipatov relation connecting the 
fragmentation functions and the distribution functions PUI, EDI ET|, we have 



D^(z) oc zq-{z), 



(34) 



where qTt{z) is the quark distribution of finding a quark q with momentum fraction 
X = z in the pion. We also know from the pQCD based analysis |^l], ^ that the 
X — 1 behavior of the quark distribution in the hadron h satisfies the counting rule 



q\x) oc (l-x)2"=-\ 



(35) 
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where rig means the spectator quark number. For favored quark fragmentation we 
notice that 

D{z)oc{l-z) (36) 
because there is only one spectator quark, whereas it is 

D{z)(x{l-z)^ (37) 

for unfavored quark fragmentation because the spectator quark number rig = 3 in this 
case. Thus we may assume 

D{z) = (1 - zyD{z) (38) 

in this extreme situation. To reflect the possible uncertainties in the parametrization 
of fragmentation functions, we also assume three options for the relation between 
favored and unfavored fragmentation functions 

D{z) = {l-zrD{z), (39) 



with = 2, 3, 4, respectively, to compare with the new parametrization ( p3|) described 
above. 

The so called Collins fragmentation function H^^^\z), which describes the tran- 
sition of a transversely polarized quark into a pion, has not been systematically mea- 
sured yet, and it is also theoretically not well known. However, there has been a so 
called Collins parametrization |^ of this fragmentation function, and we will adopt 
this parametrization in this paper. It is suggested that 

,40) 



with Mc being a typical hadronic scale around 0.3 1 GeV. Assuming a Gaussian 
type of the quark transverse momentum dependence in the unpolarized fragmentation 
function 



one obtains 



Dt{z,z'kl) = Dt{z)^exp{-R'kl), (41) 



.,,J|(._M^.V;d.^J, ,42) 
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where = j (-P^_|_), and (-P^_|_) = {kjj is the mean-square momentum the hadron 
acquires in the quark fragmentation process. We will set the parameters Mc = 0.7 
GeV and {Ph±) = (0.44)^ GeV^ as they are consistent with the spin asymmetry- 
measured at HERMES [H. 

In fact, this fragmentation function controls the z- and P/i^-dependence of the 
azimuthal asymmetries, therefore the z- and P/i^-dependence will not be considered 
in this paper. We will mainly aim at the physics of the transversity distributions, 
which controls the x-dependence of the azimuthal asymmetries. 



3 Numerical Calculations 

With the kinematics of semi-inclusive process and the formulae for the azimuthal 
spin asymmetries of pion production, and also with the quark distribution functions 
and fragmentation functions given in the above section, we can perform numerical 
calculations under different assumptions and options. We will compare the numerical 
calculations of the azimuthal asymmetries with the HERMES experimental data [||, § 
constrained by the following experimental cuts 

1 < < 15 GeV^ W > 2GeV, 0.2 < y < 0.85, 0.2 < z < 0.7. (43) 

For the and W used in the integration over y and z, we use the relations 

Q^ = sxy, = sy{l-x)+M^, (44) 

with s = 2M E = 51.8 GeV^ in the HERMES experiments. For the unpolarized quark 
distributions, we use the CTEQ parametrization (CTEQ5 set 1) as input for both 



the quark diquark model and pQCD model quark distribution functions described in 
the above section. 

3.1 Different approaches of distribution functions and frag- 
mentation functions 

We first check the difference of the azimuthal asymmetries between different ap- 
proaches of distribution functions and fragmentation functions. In this situation we 
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take into account the contribution from both the favored and unfavored fragmenta- 
tion functions, i.e., with the new parametrization of both D{z) and D{z) in Eq. (^) 
as inputs. The numerical results for Leading Approach, Approach 1, and Approach 2 
are given in Fig. |^. We find that both Approach 1 and Approach 2 differ significantly 
from Leading Approach, and this implies that contribution from the distribution func- 
tions hhix) and hi[^\x) are not negligible in the available HERMES experimental 
region. We also found that the magnitude of the azimuthal asymmetries from Leading 
Approach is below the data, thus Leading Approach cannot describe the available Tr~^ 
and 77° data. However, both Approach 1 and Approach 2 with the two model transver- 
sity distributions are consistent with the available data, except for the data point at 
X = 0.26 for TT^ production. But the predictions from Approach 1 and Approach 2 
have a large difference at large x, and we expect that further experimental data at 
large x will be able to distinguish between the two approaches. We can also arrive at 
a following conclusion: although the HERMES azimuthal asymmetries for longitudi- 
nally polarized target can provide a rough estimate of the quark transversity 
distributions by comparing theoretical calculations with the experimental data, they 
do not make a direct measurement of the quark transversity distributions since the 
contribution from other quark distribution functions might be also sizable. 

3.2 Different cases of favored and unfavored fragmentation 

We next check the contribution from the unfavored fragmentation functions D{z). 
For this purpose we compare the calculations in two cases: with both favored and 
unfavored fragmentation functions as Case 1, and with only favored fragmentation 
functions as Case 2. In order to simplify the discussion, we will only perform calcula- 
tions for Approach 2, and the results are given in Fig. |^. We find that there is almost 
no difference between the calculated results in the two cases for vr"*" and 7r° produc- 
tion, and this supports the favored fragmentation dominance assumption for the vr'^ 
and tt'' production, as was adopted in the literature. However, the situation is quite 
different in the two cases for tt^ production. This can be seen from Fig. ^ (c) and (f), 
where the calculated results for the two cases differ significantly. This implies that 
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Figure 2: The azimuthal asymmetries A^'^l^'^ for semi-inclusive vr"*", vr", and vr" produc- 
tion in deep inelastic scattering of unpolarized positron on the longitudinally polarized 
proton target. The upper row corresponds to (a) tt^, (b) 7r°, and (c) tt^ in the quark 
diquark model, and the lower row corresponds to (d) tt~^, (e) vr*^, and (f) tt~ in the 
pQCD based analysis. The dotted, dashed, and solid curves correspond to the cal- 
culated results for Leading Approach, Approach 1, and Approach 2, respectively, of 
different assumptions of distribution functions and fragmentation functions. The data 
are taken from the HERMES experiments [jl|, 0. 



the unfavored fragmentation functions have a large contribution to tt~ production. 
This can be easily understood, since the most important unfavored process for tt^ 
production is u ^ . The big difference between Case 1 and Case 2 shows that the 
u ^ process is not negligible due to the w-dominance of the quark distributions 
in the proton target at large x. 

It is interesting to notice that, although there is clear evidence for the non-zero 
azimuthal asymmetry for 7r+ and 7r° production, the available data for vr^ produc- 
tion are consistent with no azimuthal asymmetry. The data are consistent with the 



calculated results with only favored fragmentation functions, i.e.. Case 2 ^ 
However, the calculated results of Case 1 predict non-zero azimuthal asymmetry for 
TT^ production in both of the quark diquark model and the pQCD based analysis. 
From Fig. ^ (c) and (f) we find that the data point for vr" production at x = 0.26 
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Figure 3: The azimuthal asymmetries A^'^l^'^ for semi-inclusive vr"*", vr", and vr" produc- 
tion in deep inelastic scattering of unpolarized positron on the longitudinally polarized 
proton target. The upper row corresponds to (a) vr^, (b) 7r°, and (c) in the quark 
diquark model, and the lower row corresponds to (d) vr"*", (e) 7r°, and (f) tt~ in the 
pQCD based analysis. The solid and dashed curves correspond to the calculated re- 
sults for Case 1 with both favored and unfavored fragmentation and Case 2 with only 
favored fragmentation, respectively. The curves of the two cases overlap with each 
other for vr^ production. 



seems to be in disagreement with this prediction. We suggest more analysis on vr 
production at the x > 0.2 region in order to test this prediction. 

We need to examine the possible ambiguities from the uncertainties of the un- 
favored fragmentation functions. For this purpose we choose the three options of 
D{z) = (1 — z)"'D{z) with n = 2,3,4, respectively, and present the calculated results 
in Fig. ^ with cuts 0.2 < z < 0.7, and in Fig. |] with cuts 0.6 < z < 0.8, respectively. 
For the available data with cuts 0.2 < z < 0.7, the three different options of the unfa- 
vored fragmentation functions do not have big difference, even for tt^ production, as 
shown in Fig. However, for the cuts with higher z, where the unfavored fragmenta- 
tion is more suppressed than the favored ones for larger n, there is a large difference 
between the calculated results for vr" production with the three options, as shown in 
Fig. |[ Also there is a big difference between the predictions for vr" production with 
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the two different models of transversity distributions. Tlierefore tlie vr^ production 
at large x and large z is sensitive to the uncertainties of the unfavored fragmentation 
functions, and this is a region deserving careful experimental studies. 
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Figure 4: The azimuthal asymmetries A^^ for semi-inclusive vr"*", vr", and 7r~ produc- 
tion in deep inelastic scattering of unpolarized positron on the longitudinally polarized 
proton target. The upper row corresponds to (a) vr"*", (b) 7r°, and (c) 7r~ in the quark 
diquark model, and the lower row corresponds to (d) 7r+, (e) vr", and (f) 7i~ in the 
pQCD based analysis. The solid, dashed, and dotted curves correspond to the cal- 
culated results for three different options of unfavored fragmentation functions with 
n = 2,3, 4, respectively. The cuts for z are 0.2 < z < 0.7. 



With the contribution from only the favored fragmentation functions considered, 
it is found ||20| that there is a big difference between the azimuthal asymmetries 
for TT^ production predicted by the quark diquark model and by the pQCD based 
analysis, as can be seen from the dashed curves in Fig. ^ (c) and (f). However, after 
taking into account the unfavored fragmentation processes, we find a much reduced 
difference between the two different model predictions, see the solid curves in Fig. ^ 
(c) and (f). Thus a clear distinction between the two different predictions will require 
very high precision measurement in the region with both large x and large z for the 
proton target, and this is difficulty to access experimentally. It also reminds us that 
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Figure 5: The same as Fig. ^, but with cuts 0.6 < z < 0.8. 

the 7r+, 7r°, and vr" production from the proton target are mainly controlled by the u 
quark fragmentation, due to the u quark dominance in the proton target at large x. 
Thus the measured azimuthal asymmetries by the HERMES collaboration [|T], 0] are 
mainly controlled by the u quark transversity distribution 5u{x), and they are not 
sensitive to the d quark transversity distribution 5d{x). 

3.3 Deuteron and neutron targets 

We now consider the possibility to distinguish between different models by using the 
deuteron as target, where there are equal numbers of u and d quarks from isopin 
symmetry. To calculate the azimuthal asymmetries for the deuteron target, we still 
use the formalism in Sec. II, with the target polarization Sd = 0.75 multiplied by an 
correction factor 

(l - I^d) , (45) 

where ud = 0.05 ± 0.01 is the probability of the deuteron to be in the D-state. 
The azimutahl asymmetries for vr"^, 7r°, and tt~ production are calculated for differ- 
ent approaches of distribution functions and fragmentation functions, and with the 
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quark transversity distributions of the quark diquark model and the pQCD based 
analysis, respectively. For the unpolarized fragmentation functions, we use the new 
parametrization Eq. (^3]), and for the kinematical cuts we use the same as those in 
the HERMES experiments. The results are given in Fig. ^. From there we find 
that the deuteron is also not very sensitive to the two different sets of transversity 
distributions. 




0.5 1 0.5 ' 0.5 

XXX 



Figure 6: The azimuthal asymmetries A^l for semi-inclusive vr^, vr*^, and vr^ produc- 
tion in deep inelastic scattering of unpolarized positron on the longitudinally polarized 
deuteron target. The upper row corresponds to (a) 7r+, (b) 7r°, and (c) vr" in the quark 
diquark model, and the lower row corresponds to (d) vr"*", (e) vr*^, and (f) vr" in the 
pQCD based analysis. The dotted, dashed, and solid curves correspond to the cal- 
culated results for Leading Approach, Approach 1, and Approach 2, respectively, of 
different assumptions of quark distributions and fragmentation functions. 



We now check whether the situation can be improved if we use the neutron as 
target. This can be experimentally accessed by measuring the numerators and denom- 
inators of the spin asymmetries (|10D from proton and deuteron targets respectively, 
and then extract the numerator and denominator for the neutron contribution to cal- 
culate the spin symmetry of the neutron. The calculated results are given in Fig. ^ 
We find that there is a somewhat bigger difference between the predictions of the 
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Figure 7: The same as Fig. but with the neutron target. 

quark diquark model and the pQCD based analysis. In order to understand this 
behavior, we consider the vr"^ production for a neutron target, shown in Fig. |^ (a) 
and (d), and compare them with the vr" production for a proton target, shown in 
Fig. 1^ (c) and (f). For the proton case, the favored fragmentation process is d ^ tt~ , 
while the unfavored process is -u ^ vr^. The unfavored process is important due to 
the dominance of u quark distribution over d inside the proton target at large x. For 
the neutron case, the favored fragmentation process is u —>■ vr"^, while the unfavored 
process is d ^ tt^ . The unfavored process is suppressed by a factor of 4 (from the 
squared charge) compared to the corresponding process of the proton target, while 
the favored process is enhanced by a factor of 4 compared to the corresponding pro- 
cess of the proton target. Thus the favored process u relative to the unfavored 
process d ^ for the neutron target is enhanced by a factor of 16 compared to the 
corresponding favored process d ^ relative to the unfavored process u ^ tt^ for 
the proton. This implied that the favored process of 7r+ production for the neutron 
target plays a more important role than the corresponding process of tt~ production 
for the proton target, where the unfavored process is more important and cannot be 
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neglected. Thus we suggest to distinguish between the quark diquark model and the 
pQCD based analysis by measuring 7r+ production from a neutron target. The vr"*" 
production for a neutron target is sensitive to the u quark transversity distribution of 
the neutron, which is essentially the d quark transversity distribution of the proton. 
Thus combination of pion production from both the proton and neutron targets can 
measure the flavor structure of the quark transversity distributions. 



3.4 The azimuthal asymmetries of 

For the azimuthal asymmetries with weighting function sin 20, i.e., A^'^J^'^, there is 
only one quark distribution function hfl^\x) involved in the numerator of (|ll|)- 
both the Leading Approach and Approach 1, this distribution function is assumed 
to be zero. Therefore the azimuthal asymmetries A^^J^'^ can only be calculated in 
Approach 2, and we present the results for the proton, deuteron, and neutron targets 
in Figs. ^ and 0, respectively, in the two different models of quark transversity 
distributions. We find that the magnitude of these asymmetries are not so small in 
the two models, thus it is possible to measure them with high precision measurements. 
In the case of a proton target, the calculated results are similar for the two models 
with the available HERMES kinematical cuts, and they are both compatible with the 
available HERMES data. The calculated results are also similar in the two models 
in the case of the deteron target. However, larger differences between the predictions 
of the two models is found in the case of vr"*" production from a neutron target. Thus 
precision measurement of A^^J^'^ from nucleon targets may provide us with information 
on the distribution function hf[^\x) and its flavor structure. 
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Figure 8: The azimuthal asymmetries for semi-inclusive tt"*", 7r°, and 7r~ pro- 

duction in deep inelastic scattering of unpolarized positron on the longitudinally 
polarized proton target (with polarization S = 0.86). The upper row corresponds 
to (a) TT"^, (b) 7r°, and (c) in the quark diquark model, and the lower row cor- 
responds to (d) 7r+, (e) 7r°, and (f) 7i~ in the pQCD based analysis. The solid and 
dashed curves correspond to the calculated results for Case 1 with both favored and 
unfavored fragmentation and Case 2 with only favored fragmentation, respectively. 

4 Prediction of Pion Azimuthal Asymmetries for 
Transversely Polarized Target 

From the above discussions, we find that although the azimuthal asymmetries off lon- 
gitudinally polarized target may serve as a measurement of the transversity distribu- 
tions, there are still some contribution from several unmeasured distribution functions 
and fragmentation functions. Their contribution are not negligible in the HERMES 
kinematical region, and also there are some ambiguities concerning the size of their 
contribution from different assumptions, as we can see from the different predictions 
between Approach 1 and Approach 2 in the last section. Therefore it is necessary 
to look for more clean processes for a direct measurement of the quark transversity. 
In fact, it has been known that the azimuthal asymmetries off transversely polarized 
target are directly connected to the Leading Approach term we discussed for the 
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Figure 9: The same as Fig. ^ but for the deuteron target with polarization S = 0.75. 



longitudinally polarized situation. There has been some preliminary results from the 
SMC experiment |^3[|, indicating evidence of non-zero azimuthal asymmetry for pion 
production from a transversely polarized target. Also the HERMES collaboration is 
planning to measure the azimuthal asymmetries of pion production from transversely 
polarized targets in the near future . It is thus necessary to make clear predictions 
on the quantities they will measure by using the two sets of transversity distributions 
given by the quark diquark model and the pQCD based analysis. 

More specifically, for an unpolarized lepton beam and a transversely polarized 
target, the following weighted asymmetry provides access to the quark transversity 



distribution 10 



ATix,y,z) 



/d<^7d2^/.^^sin 



d(T^ — da-^ 



(46) 



where t (J.) denotes target up (down) transverse polarization. The azimuthal angles 
are defined in the transverse space giving the orientation of the lepton plane (0') and 
the orientation of the hadron plane (0'^ = (t>h ~ 00 oi spin vector {(f)\ = (pg — </)'■) 
with respect to the lepton plane. The asymmetry (EB) can be calculated from the 
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Figure 10: The same as Fig. but for the neutron target with polarization S = 0.75. 



distribution functions and fragmentation functions by ||T0| 

ATix,y,z) = St 



(47) 



(l-y + y2/2)j:^e^^q{x)Dl{zy 
where St is the target polarization. Assuming St = 0.75 and the same kinemat- 
ical cuts as in the MERMES experiments, and with all of the quantities discussed 
as before, we can calculate the asymmetry At with the two different sets of quark 
transversity distributions. The integration over y and z should be performed for 
both the numerator and denominator. The predictions for a proton target are given 
in Fig.|lT], from where we notice that the azimuthal asymmetries of vr"^, n^, and tt~ 
production from a proton target are mainly controlled by the u quark transversity 
distribution, and they can serve to measure 6u{x). The predictions for pion produc- 
tion from a deuteron target, given in Fig. |12], are also found to have small differences 
between the two different models. However, the azimuthal asymmetries from a neu- 
tron target are predicted to have a large difference between the quark diquark model 



and the pQCD based analysis, as shown in Fig. |T3|. Thus they can serve to measure 
the flavor structure of quark transversity distributions, i.e., one can measure both 
6u{x) and 6d{x) by using the azimuthal asymmetries from both proton and neutron 
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Figure 11: The azimuthal asymmetries {At{x)) for semi-inclusive tt"*", 7r°, and tt~ 
production in deep inelastic scattering of unpolarized positron on the transversely 
polarized proton target. The upper row corresponds to (a) 7r+, (b) tt", and (c) 7r~ 
in the quark diquark model, and the lower row corresponds to (d) tt"*", (e) 7r°, and 
(f) 7r~ in the pQCD based analysis. The solid and dashed curves correspond to the 
calculated results for Case 1 with both favored and unfavored fragmentation and Case 
2 with only favored fragmentation, respectively. 



5 Summary and Conclusions 

In this paper, we analyzed in detail the azimuthal spin asymmetries of pion produc- 
tion, in semi-inclusive deep-inelastic scattering of unpolarized charged lepton beams 
on longitudinally and transversely polarized nucleon targets. Various assumptions 

and approximations of the distribution functions and fragmentation functions were 
examined. It was found that different approaches of distribution functions and frag- 
mentation functions lead to different predictions of the azimuthal asymmetries in the 
available HERMES kinematical region. This means that the effects from unknown 
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Figure 12: The same as Fig. [T^, but for the deuteron target. 

distribution functions may have sizable effects on the azimuthal asymmetries of the 
HERMES experiments, thus one needs to consider them before attributing the HER- 
MES data as a measurement of transversity distributions. We also found that the 
predictions of tt~ production from the proton target are quite different for two cases 
of both favored and unfavored fragmentation, and of only favored fragmentation. 
This means that the unfavored fragmentation functions play an important role for 
TT^ production from a proton target, due to the dominance of u quarks inside the 
proton target, thus the unfavored u tt~ fragmentation are not negligible. This 
point was not considered in the available studies of azimuthal asymmetries. The pion 
production from the proton target is most suitable to study the u quark transver- 
sity distribution. In combination with the pion production from a neutron target, 
it is possible to measure both of the u and d quark transversity distributions. Thus 
we suggest to measure the azimuthal symmetries of pion production from the neu- 
tron target, in order to test different model predictions. Predictions of azimuthal 
asymmetries of pion production in semi-inclusive deep-inelastic scattering of unpolar- 
ized charged lepton beams on transversely polarized nucleon targets are also given. 
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Figure 13: The same as Fig. |TT], but for the neutron target. 

and different predictions between two different models are given for pion production 
from a neutron target. This study will be useful for designing experiments aiming at 
measuring transversity distributions through pion production in semi-inclusive deep- 
inelastic scattering of unpolarized beams on longitudinally and transversely polarized 
targets. It is also useful for theoretical studies aiming at extracting transversity distri- 
butions from experimental observation of azimuthal asymmetries of pion production 
in semi-inclusive deep-inelastic scattering. 
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